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Abstract 

Using the ellipsoidal model for the density configuration, we calculate the 
equilibrium sequence of the corotating binary stars of the polytropic equa- 
tion of state in the first post-Newtonian approximation of general relativity. 
After we calibrate this model by comparing with previous numerical results, 
we perform the stability analysis by calculating the energy and the angular 
momentum of the system as a function of the orbital separation. We find 
that the orbital angular velocity at the energy and/or momentum minimum 
increases with the increase of the compactness of each star, and this fact holds 
irrespective of the polytropic index. These features agree with those in previ- 
ous numerical works. We also show that due to the influence of the tidal field 
from the companion star, the central density of each star slightly decreases. 
PACS number (s): 04.30.Db, 04.25.Nx, 04.40.Dg 
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I. INTRODUCTION 



In a previous paper |l| , we analytically obtain the equilibrium configuration of corotating 
binary stars(i.e., the solution of the Darwin problem) of incompressible fluid in the first 
post-Newtonian(PN) approximation. We calculate the energy and angular momentum as a 
function of the orbital separation, and investigate the location of their minimums(we call the 
orbit at the minimum the innermost stable corotating circular orbit (ISCCO)Q). It is found 
that the angular velocity at the ISCCO increases with the increase of the compactness of 
the star; i.e., the general relativistic(GR) effect makes the angular velocity at the ISCCO 
larger. 

In this paper, we try to extend the previous work to the case of the compressible fluid star 
using the ellipsoidal model for the density profile(the so called ellipsoidal approximation) 
proposed by Lai, Rasio and Shapiro ||. The ellipsoidal approximation may be used for 
stars of the sufficiently stiff equation of state(EOS) of a large adiabatic index. Thus, this 
method may make it possible to investigate the ISCCO for realistic binary neutron stars 
which are made of the stiff EOS of the large adiabatic index r ~ 2-3 |. In previous 
papers, one of us(M.S.) || numerically calculates the equilibrium state of binary neutron 
stars of the polytropic EOS with T = 2 and 3 in the PN approximation. We will show 
that by a simple semi-analytic calculation, the results obtained in previous papers can be 
qualitatively explained well. Also, we investigate the influence due to the GR tidal effect in 
binary neutron stars. 

In section II, we show the basic equations to derive the angular velocity for the corotating 
binary stars. In section III, we introduce the ellipsoidal approximation assuming the poly- 
tropic EOS. Using the ellipsoidal approximation, in section IV, we calculate the equations 



*The ISCCO is the orbit where the secular instability sets in, and different from the innermost 
stable circular orbit(ISCO) where the dynamical instability sets in. The ISCO will locate inside 
the ISCCO §. 
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for the energy and angular momentum. In section V, we derive the scalar virial relation in 
the PN approximation, which is used to determine the size of the star. In section VI, we 
show equilibrium sequences of binary stars for various polytropic indices. Paying attention 
to the energy and angular momentum as a function of the orbital separation, we investigate 
the ISCCO. We also show that due to the influence of the PN tidal field from the companion 
star, the central density of each star slightly decreases. Section VII is devoted to summary. 

Throughout this paper, we take the units of G = 1, where G is the gravitational constant, 
and use c as the light velocity. Latin indices i, j, k, I, ■■ take 1 to 3, and 8^ denotes the 
kronecker's delta. M denotes the Newtonian mass for each star of the binary as 

M = J pd 3 x. (1.1) 

The quadrupole moment and its tracefree part of each star are defined as 

hi = J px\d z x, and l u = I u - - ^ I u . (1.2) 



II. FORMULATION 

In this section, we show the basic equations deriving the equations of the angular velocity 
and the mass for corotating binary stars in the first PN approximation with no restriction to 
a special model because these quantities can be written in a simple form without specifying 
a density configuration. To construct the equilibrium sequence and argue the stability, we 
additionally need the equations for the energy, the angular momentum, and the scalar virial 
relation in the PN approximation. We derive them after we specify the model for the density 
configuration because they are not written in a simple form. 

We consider corotating binary stars of equal Newtonian masses(M! = M 2 = M), whose 
coordinate separation is R, in the first PN approximation. Here, the gauge condition is 
the standard PN one ||. We adopt the coordinates where the center mass of a star (star 
1) locates at the origin and the other one(star 2) locates at (xi,x 2 ,x 3 ) = (—R,0,0). As 
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for the density configuration, in this section, we only assume that the star has the triplane 
symmetry with respect to the symmetric plane(i.e., for star 1, p(x,y,z) = p(—x,y,z) = 
p(x, —y, z) = p(x, y, —z)). In the following, when the integral symbol appears, it means that 
we perform the integral inside the star 1. 

In deriving the equations for various quantities, we assume R> a > M/c 2 , where a is 
the typical size of the star: We calculate the terms up to O((ao/-R) 3 ) relative to the lowest 
order for the angular velocity and the angular momentum, and the terms to 0((ciq/ R) 4 ) for 
the energy. Then, the effects by the quadrupole moment of each star to the equilibrium 
state are taken into account consistently. 

We adopt the polytropic EOS as 

p = Kp T = (r - l)pe, T = 1 + -, (2.1) 

n 

where P, p, e, n and K are the pressure, the mass density, the specific internal energy, the 
polytropic index and the polytropic constant, respectively. Then, in the first PN approxi- 
mation, the integrated form of the Euler equation is written as |J [f| 

1 / ^ 2 



= U ~ ff + + ^( 2uj2f/ - X ° + ^)} fi2 + + constant, (2.2) 

where w 2 = (x x + R/2) 2 + x\ and K' = K(n + 1). 

The angular velocity is obtained from the first tensor virial(TV) relation as 

= J (l - = J pU^x + *MQ 2 + I [- J P X ^x 

+n 2 x J p(2tu 2 U A + A Xl U + 2RU - X nA + ^ A )d 3 x 



i MR 



3\ -\ 



+ l4(-(3/ 11 + I 22 , 



(2.3) 



where , k denotes the partial derivative with respect to x^. denotes the angular velocity 
in the Newtonian order solution as HI 



To determine the orbital angular velocity up to the PN order, we need to perform the 
following integrals: 

I l = J p Ud 3 x, (2.5) 

J P U A d 3 x, (2.6) 

J pU Xl d 3 x, (2.7) 

J pU lXl d 3 x, (2.8) 

J pU lX 2 k d 3 x, (2.9) 

J P X 0A d 3 x, (2.10) 

J P X nA d 3 x, (2.11) 

J pP v ,id 3 x. (2.12) 



h 
h 
h 
h 
h 
h 



In the following subsections, we separately show the equations for the gravitational potentials 
and their solutions to evaluate I x ~ I 8 . 

A. Equation and solution for U and integrals including U: l\ ~ 1$ 

Equation for the Newtonian potential U is 

AU = -Anp. (2.13) 

Here, there are two kinds of contribution for U at star 1; one is the contribution from star 
1 itself and the other is that from star 2. Hence, we write U as 

U = U 1 ~* 1 + U 2 ^\ (2.14) 

where U l ^i denotes the potential generated by star i at star j. Here, U 2 ^ 1 is generally 
expressed as 

TT 2^i M( x 1 2x\-x\-x 2 -2x 3 + 3 Xl (x 2 + x 2 ) ^ D _ 4 ^ 
U = RV-R + 2^ + 2^ + ° {R ] ) 



+ 
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On the other hand, we cannot specify U 1 ^ 1 without the density configuration. Thus, we 
define the following integrals M; 



W = ^J P U 1 ^ 1 d 3 x, 
Ma = I pU^Xid^x, 



(2.16) 
(2.17) 



where we have the relation J2i-M.u = —W ||. Using the solution for U, we get 

M 2 9MJ U 6 
M 

h = -— 2 In + 0{B7±), 

h = M 11 + '^I 11 + 0(R~ 5 ), 

h = -^- 2 I kk + 0(R- 4 ). 

Note that fl^ is derived using I 2 . 



(2.18) 
(2.19) 
(2.20) 
(2.21) 
(2.22) 



B. Equation and solution for Xq and integrals including Xq: Iq 



As in the case of U, X is divided into two parts as X Q = X^ 1 + Xq^ 1 . We also write Iq 
as Iq^ 1 + where I}? 1 and If - * 1 correspond to the contributions from Xq^ 1 and X^ 1 , 



respectively. Equations for Xq^ 1 and X^ 1 are 



47rp 



AXq 2 ^ 1 = top 



2C/ 1 ^ 1 + 2U 2 ^ 1 +e + 



2U 2^2 + 2U i-*2 + £ + 



3P 
P 

3_P 

P 



[2.23) 
[2.2A) 



To evaluate Jg , we need only the odd function of x% in pj? , and the terms required are 



p^i = _ 2p ^i ( M + M 2Xl 3X \ ^ + 



r + 0(i?- 4 )). 



(2.25) 



Then, the terms of Xq required for 1$ are 

X^ 1 = + ^) A + f(2^m - 3£ 122 - 3D 133 ) + 0(iT 6 ) 



(2.26) 



where Di and D\ kk are the solutions of 



ADi = -Airpxi, (2.27) 
AD lkk = -Anpxixl (2.28) 



Using the relations, 



J pD iA d 3 x = -Ma, (2.29) 
J pD lkkA d 3 x = - J pU^ 1 x 1 x 2 k d 3 x = -M ukk , (2.30) 

Jg^ 1 is evaluated to be 



1(3 = " ll? + "R^J^ 11 " i^ 2 -^ 1111 " 3jMl122 " ^i^) + 0(R- b ). (2.31) 
On the other hand, the solution of Xq^ 1 is written as 

xr = f (i - 1 + + =M±j|^±j£l£i + (^)) 

+ l( 1 -^ + °( fl - 2 )) + W i ( 1 -^ + °M' < 232 » 

where 

M x = - J p^d 3 x = -AW - (3r - 2)Ui " 2 + + 0(7T 5 ), (2.33) 

= - / = + 0(iT 4 ), (2.34) 
J-xn = / Px^M -x\- x\)d 3 x 

= -2+ uu - (3r - 2)^ e n - ^J^n + 0(iT 3 ), (2.35) 



and 



[/; = y perf 3 a;, (2.36) 

rf i/ii = , Iukk = J pU x k d x, (2.37) 

T _ 2/ e ll — / £ 22 — / £ 33 T 

+ell — o , Jefefe 



J pex 2 k d 3 x. (2.38) 



Thus, /g 2- * 1 is evaluated to be 
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I^ 1 = -M, 
6 v R 2 



M x , 2d x , 9M x ^ n , 9^ X11 



i? 3 2Mi? 4 



+ 



2i? 4 



= M 



4W + (3r -2)Ui 2M 2 2AM 4-u AMI U 



R 2 



R 3 



R 5 



R 5 



+ 2^{( AW + (3r " 2)f/i )l7 + 2 + uu + (3r " 2) ^ n } + 0(jR_6) 



(2.39) 



C. Equation and solution for Xq and integral including Xq: I? 



Equation for X n is 



AX n = 8np(xl + xj + Rx! + 



4 



(2.40) 



The solution is also written as X^ 1 + Xq* 1 , where 

X^f 1 = -2 (Dn + D 22 + RD 1 + ^V^ 1 ) , (2.41) 
XI- 1 = -^U^ + 2|(l - |l + 0(iT 2 )) - 2(711 + /22) (l - I + 0(R^)) , (2.42) 

and _D fcfc is the solution for AD kk = —\^px\. Using these solutions, we soon get I 7 as 

h = 2RMu + y(M + - —(In - I 22 ) + 0(i?- 3 )J. (2.43) 



D. Equation and solution for (3^ and integral including (3^: Is 



Definition of P v is 

P v = ~(x 1 P 1 + x 2 P 2 + ^P 1 > j 



where 



R^ 2 
2 



AP 1 = -4^(^ + 1), 
AP 2 = -4vrpx 2 . 



^2.44) 



(2.45) 
(2.46) 
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Pi is also written as -P/^ 1 + -Pj 2_>1 , where 



P^ 1 = D x + *U^\ (2.47) 
P^ 1 = D 2 , (2.48) 



and 



P^ = *f 2 (l - ^ + 2 3 | + 0(i?- 3 )) - f (2.49) 

PT 1 -^-^ + 0(RT*)\ (2.50) 



i? 2 V R R 2 



fl 2 ^ 1 can be written to the explicit form as 



2 ^ MR/ 7m 14x 2 + 9^ + 7^ x 1 (14x 2 + 13^ + 7 3 :|) n ,„_ 4 ^ 
^ = ~1 7+ 1T 2R 2 + 2^ + ° {R ] ) 

Then, after straightforward calculations, J 8 is evaluated as 

/» <j 1/2 1 / 

/ 8 = _-(7W - M 22 ) + — + ^( 14/ n + 9/ 22 + 7/ 33 ) + 0(iT 3 ). (2.52) 



E. ft 2 



Gathering the results obtained in previous subsections, f2 2 is calculated as 



Q 2 



2M 



1 r 9 

1 + M^r w+ (3r " 2)f/i + Maa ~ Mu ~ ~ 

(28/n - I6/22 - 5/33) + 0(iT 4 )} 



+ 



i? 5 



2i? 3 

1 / sr — 2 

1 + ^ (w - A^n - M22 + -^—u, + o(/r 2 ; 



+ ^e( 18 +uu + 9(3r - 2)^11 - 4M im + QM U22 + 6Mn 33 ). (2.53) 

To transform this equation into a more simple form, we use the scalar, the second and the 
fourth tensor virial equations in the Newtonian order as follows; 
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3(r - i)Ui = w 



3M 



R 3 



/n-54(/n + / 22 ) + 0(iT 5 ), 



Mll = _ (r _ - nlh 



2MI 



11 



R 3 



0(R 



-5\ 



M 22 = -(r - l)Z7i - tt 2 / 22 + ^ + 0(R~ 5 ) 

x 33 = -(r - m + — ^ + o(R~ 5 ), 



133 — 
2 M mi 



3.M 



1122 



i? 3 

3^1133 



-9(T-l)+ £ll + 0(R 



-3\ 



(2.54) 
(2.55) 
(2.56) 
(2.57) 
(2.58) 



Then, 



2M 



R 3 
18Ju 
i? 5 



1 r 9M 2 



Mc 2 



M 



47? 2i? 3 



(28/n - 14/22 - 9/33) + o(r- 4 : 



1 + ?( h ?f +0 ^ 2 : 



+ 



c 2 # 5 



18^t;ii + 9(5r-4)^ 



Ell 



^2.59) 



We note that in the incompressible case(p =constant and T — > 00), the configuration of each 
star is the ellipsoidal figure of its axial length ai, o 2 , and 03, and we have relations 

2MnpA 



W 



Ui = 0, 



where 



A = aia 2 a 3 / 
Jo 



2 

a x a 2 a3 



16npA 
~~ — 21 — 

F-i)^ = ^, 



^ (of + u) (a 2 + uj (a 2 + w) 
eft 



(2.60) 
(2.61) 
(2.62) 
(2.63) 



= alA (a 2 ,a 3 ) 



(2.64) 



and a 2 = a 2 /ai and a 3 = 0.3/01. Then, Q 2 reduces to that derived in a previous paper 



F. Mass and center of mass 

The conserved mass is defined as [|J 
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M * = J d 3 xp* = J d 3 xp 



1 + 



1 (v 



c 2 \2 



+ 3U 



M 



1 + 



1 /6W 13M 



+ ~nr + 77^( 45J1 ii + 4/ s) + 0(R~ 



c 2 V M AR AR 3 
where I s = In + 122- The PPN mass || is slightly different from M* as 

1 (v 2 „_ f/ 1 ^ 1 



(2.65) 



Mi 



PPN 



.1/ 



1 f5W + U { 13M 1 _ r . ^. 

+ 7^(45^ 11 + 8/ s ) + 0(i? 



(2.66) 



M 4R 4.R 3 

The advantage of the PPN mass is that if we use it, Q 2 can be denoted only by Mppn for the 
limit hk/R 3 —> |§. However, when I^k is large enough, it is not the conserved quantity 
any more, and not so useful. 

We define the center of mass by the conserved mass density as 



x- 



1 

M* 



^2.67) 



and the x\ coordinate of the center of mass for star 1 becomes 



1h 



r + 0(R- 



R 2 

Note that in the PPN formalism ||, the center of mass of each star is defined as 
1 



(2.68) 



x 



PPN 



Mi 



PPN 



d xpx 1 



1 (v 2 U 1 ^ 1 O 2 

1 + Ay +£ + 3U - — + 2 {xl+xl 



(2.69) 



and the results are the same. From Eq. (|2.68 ) we define the orbital separation as 



R* — R 



PPN 



R<1 + 



4/n 
' R 3 



+ 0(R~ 



(2.70) 



III. THE ELLIPSOIDAL APPROXIMATION 

In this section, we introduce the ellipsoidal approximation for the density configuration 
proposed by Lai, Rasio and Shapiro(LRS) 0. In this approximation, the configuration of 
each star is assumed to be the ellipsoidal figure of its axial length oj, ci2 and 03. Then, the 
density profile is assumed to be 

11 



where p p (£) has the density profile of the spherical polytrope, and p = for £ > £1. Because 
of the polytropic relation, P and £ are also assumed to depend only on £. 

Here, there is no solid basis to justify this model: i.e., in this approximation, we do not 
rigidly treat (a) the change of the density profile due to the PN gravity which occurs even 
for the spherical star, and (b)the deformation of the star due to the tidal gravity by the 
companion star. These facts mean that several quantities of the system calculated below 
such as the energy, the angular momentum, the central density of the star, and so on, deviate 
systematically from their exact values(see section VI). However, (1) in the Newtonian case, 
this treatment is almost exact for the incompressible fluid M |J, and (2) in a previous 
paper [T|], we find that the tidal deformation by the PN gravity is a small effect for the 
incompressible fluid star. Thus, we may expect that the ellipsoidal approximation is good 
for a sufficiently stiff EOS (small n). 

Following the standard procedure, we set 



where p c is the central density, and 9 obeys the Lane-Emden equation ||. Then, M, Ikk, 
U[, and I e kk are immediately calculated as 



(3.2) 



M = A.Tip c aia2a^ 
r _ Ma l . 



(3.3) 





(3.4) 



Ui = hKp, 



>l /n M, 



n(n + 1) 



(3.5) 



5 — n 





(3.6) 



where & = 89/ 'd£ at 



As for W and A4kk, we follow LRS and set as 




5 — n 



(3.7) 




(3.8) 
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where p = 3M/(47raia2a3), in order to guarantee that in the case of the incompressible 
ellipsoid or the spherical star, these are exact definitions. Due to the same principle, Iukk is 
written as 



lUkk 



5-nW 
3 M 



5a r , 



(3.9) 



where we make use of the fact that in the case of the spherical star and/or the incompressible 
star, U 1 ^ 1 is written as 



U 



i-»i 



5 -nW 



1 







(3.10) 



3 M\ ^\9'\ 

Thus, all quantities to calculate Q 2 and M* reduce to the available form, and they are 
rewritten as 



fi 2 



M f 2 



+ 



6Kngil 

5i? 5 



+ 



M 



2 r 



c 2 af 



r?\ 
+ 



5W + Ui 
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4/2 10i? 3 



(2? 



14a 2 . - 9a£ 



6n n q n T j,/23 — n 



5# 5 



-W 



6 



+ 7r, 



15n + 35 
6{n + 1) 



M* = M 



1 + 



M 
c 2 ai 



+ — + -^-(15gn + 4g B 
4i? 2(LR 3 



(3.11) 
(3.12) 



where gn = 2 — a 2 — cn|, q s = 1 + a|, -R = R/a>i, In = (5 — , n.)cr ri /(^i|6 l '|K n ), and 



3A 



M 2 2(5-n)a 2 a 3 ' 



-iZ N 



(3.13) 
(3.14) 



5R 3 5 

Here, to derive the expression for t/j, we use the scalar virial relation in the Newtonian order. 
In the ellipsoidal approximation, a 2 and 013, which are functions of R, are assumed to be 
determined from the tensor virial relations in the Newtonian order as M 



f - M = -k x Kpl' n M + — ^ + fi= I u , 
5 — n n /i J 

5 — n n R 3 

2iipA z al w 1, MJ33 
, 3 M = -ktKpl^M - —p. 
5 — n n R 3 



'22, 



(3.15) 
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IV. THE ENERGY AND THE ANGULAR MOMENTUM 



Next, we consider the total energy of the system. The energy for each star is calculated 
from E 1= E N + E PN /c 2 , where § § 



E N = J p (e + ^v 2 -^Uy 3 x, 
E PN = J p(jV + h -v 2 U + Tev 2 + 2eU - \\J 2 + ^ 2 )d 3 *, 



(4.1) 
(4.2) 



and v 2 = w 2 Q 2 . Even for the compressible PN Darwin problem, the expression of E-& is 
soon derived by using the quantities defined above as 



E N = U i -W + 

(M 2 
V at 
(M 2 



Q 2 


fMR 2 


r \ M l 




_3^A 


2 


\ 4 




,r~ 


R 3 J 



fi^ (W_ Kn(h\ _ 1 / 1_ K n q U 

2 V 4 5 j 2 \R BR* 



(4.3) 



2a lC 2 ~™ V 4 ' 5 

where tt 2 = Q 2 a\/M = figr + (M/cV )^pn- On the other hand, to express £<pn, we are 
required to evaluate new integrals as 



J peU^cPx, (4.4) 

J piU^fcPx, (4.5) 

J pw^U^cfx, (4.6) 

J pU^U^Wx, (4.7) 

/ pfiT^x. (4.8) 

To evaluate Ji ~ J 4 , we use the same principle in calculating ^un- Then, they become 

(4.9) 



Ji 
J* 

h 

■h 



Ji = ^Kp^Wm(^^ + n + l 

5-nW 2 f(5-n)X n 
J<i = — —I .o,Jo +n + 7 



■h 



9 Ml (l\0'\ 2 
R 2 W 5-nW 



+ 



/ s 5 



2 3 M V 



+ 1 



^ 2W '(i + w) + 5|( 5 -") H '( 5 ^fc + 1 ) +0 < ff 



(4.10) 
(4.11) 
(4.12) 
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^-if 1 ^ (4.i3) 



where 



am Jo 

On the other hand, J 5 is written as 

J 5 = ^-(M 22 - 7W )~lJ P( D ^ + te)ci 3 x 

~ J p(x 2 2 D 1A + x?D 2>2 - X!X 2 D lt2 - xix 2 D 2A )d 3 x. (4.14) 

Unfortunately, the integrand in the second and third terms cannot be expressed by 0. Hence, 
we need to introduce a new function: To evaluate Di and D 2 , we assume that the density 
profile p, which appears in the source terms for their Poisson equations, is equal to that of 
the spherical polytrope p p (£). Then, we can write D k as 

D k = Anp c a 2 4>{^)x k , (4.15) 

where a 3 = a\d 2 d^ and satisfies 

Note that the following relation holds; 

EA^^ + E^jr 1 . ( 4 - 17 ) 

k k 



Also, for the spherical star, 



2 •") // 11' 



4wi = (418) 

/,From these two equations, we have the relation 

30 + ^ = 6l^l+^ + e|. (4.19) 

Thus, <fi is solved under the boundary conditions = {i\\0'\ + l)/3 and d<fi/dt; = at £ = 0. 
Using this definition of D k , J 5 can be written as 

J, = ^ {M22 - 7W) .-JM^ M + 4)Vn , (4.20) 

15 



where 



Exchanging M 2 j '(a 1 a 2 a 3 ) 1 / 3 into (5 — n)W/3, we finally get 



d2 AW 

y (_M 22 - 7W) - -|-(5 - n)(a? + a 2 )^ 



(4.21) 



(4.22) 



Gathering the above results, we obtain the expression for Ep^ as 



E 



PN 



X n + n + 7 + 



5(5 - w) IU 2 
18 ~M 

4 V 2 2/ n 4 



2(5 -w) t/jW /^ 
3(n+l) M V 



A„ + n + 1 



+ 



1 

Z" 3 



11 



33 - lOn 25 



25 \ T 
+ y 7«J + Is 



W^~E 

9 kA Ui(fl5n 9 
4 n 



7M 3 
32i? 2 

25 5(5 -n) 20(5 - n) ^ 
y7n+ g g — 

57, 



+ 7U11 + — =7- (2(n + l)/ s + 3n^ n 
2/ n + 1 v 



+ 



M 2 



■(1O6/11 + 25/ 22 - 71/33), 



(4.23) 



16i? 4 

where A„ = (5 — n)\ n /(£ 2 \9'\ 2 ). Using the scalar and tensor virial relations in the Newtonian 
order, Ep^ is rewritten to a simple form 



E t 



PN 



5(5 - n) W 2 
18 ~M 
M ( 55 



A n + n + 7 + 



2(5 -n) UiW f~ 



+ R 

w 

+ R* 

M 2 



48i? 4 



W + -Ui 
6 2 

Jui (571 - 20 + 



7M 3 



32i? 2 



5n 
n + 1 



3(n+l) M 



25 



A„ + n + 1 



7n - — 7n + /s — 7n + 



35 



5(5 - n) 20(5 - n) rj r , 



9 



302/ n + 59/22 - 209/33 +0(R 



>-5\ 



where 
Epn ; 



/M 3 
5(5 -n 



-Ep>N (-R) , 



(4.24) 



+ - 



R 



55. 



W 2 X n + n + 7 + 



5. 



2(5 -n) 
3(n+ 1) 



6 



W + -Ui)- 



n n w r 



+ 



15i? 3 

K 



q n 5n - 20 + 



32/? 2 

5n 



-lr, 



240/? 4 



C"V"" ' n + 1 "' 2 
^y (302 + 59a 2 . - 209a 2 .) + 0(/T 5 ). 



25 U f« ,r^r: N 20(5 - n) ^ 

7n J + 5s I 357„ + 5(5 - n) 



K. 



(4.25) 
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For the incompressible case n = 0, k = 1, o"o = 2/35, 70 = 1/7, r] = 6/35, Ao = 8/35(see 
Table I), U- x = 0, and W/M = 27rpA /5, the expression of E PN , of course, agrees with that 
obtained in a previous paper 

The angular momentum for each star J\ = Jn + Jpn / c 2 is calculated from 



Jn 

JpN 



J p [v 9 (n 2 N m 2 + 6U + Te) + f3 v Q 



d 3 x, 



where v v = flzu 2 . Using J\ ~ J 5 , they are soon expressed as 



_ _ /M/i 2 
J N = I — — + Is 



Ji 



PN 



4 
i? 2 



^ f35 7 n 



• M 3/ y /2 ftjN(£), 

M 

+—(20/ n + 15/22-11/33) 
4Jx 



(4.26) 



(4.27) 



M 5 / 2 ar 1/2 fi r 



— (f>i + 5#) + -/? + 
4 v ' 2 5 



K ra g s iy f 357 n 



20 



1 _- r 2(5-n)-^(5-n)^ 
13 9 k„ 



+ 



20/? 

= M 5 / 2 a^ 1/2 fi N J PN ( J R). 



(20 + 15a 2 - 11a 2 ) 



(4.28) 



Once we determine a^(R) and a^R) from Eqs.( p.l5| ), we can calculate An, ^pn -^n, 
-Epn, Jn, and Jpn as a function of .R right now using the numerical values for £ l5 fci, 
k„, <T n , A n and shown in Table I. To evaluate the angular velocity, the energy and the 
angular momentum, however, we further need the equation for determination of oi. We will 
derive it in the next section. 



V. SCALAR VIRIAL RELATION IN THE POST-NEWTONIAN 

APPROXIMATION 

Up to this section, we have derived the equations required for determining the equilibrium 
sequence except for ai(or p c ). Here, we derive the equation for ai(or p c ) in the PN order. 
Such an equation should come from the scalar virial relation in the PN approximation, which 
is derived by 
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*v--^(*v-^ 



(5.1) 



Using the equations shown in previous sections, the left-hand side of Eq. (|5.1|) can be written 



as 



M 



n c z n + 2 



(5.2) 



In the right-hand side of Eq.( |5.1| ), many terms are soon integrated as 



R 6 cr 



rx A5 + n TjA 6M 2 \ , 

i? 2 



+^|2 Y^Mau + M ll22 ) - —W + i{I un + I U22 ) -K 2 + K 3 
+— (2/ u + 15/ 22 - 15/ 33 )] + ^^ 2 (3/n + I 22 ) + 0(iT 5 ) 



(5.3) 



but for 



#1 = E / ^irw 

i 
i 

^3 = E / p/^TW 



;r. 



(5.4) 
(5.5) 
(5.6) 



i^i is evaluated as 



Ki=-j y^ 3 xAp(x)p^ i (i/)E^^(^3 

= I J d 3 xd 3 yp(x)p 1 x - 1 (y)Y: 

i 



d 



1 



% V|x- y| 



/A^w^+E^r 1 : 



1 + - ( J l + J 'l) + 2 ( J 2 + J2 + J 4 + J4), 



(5.7) 



where we define 



i 

J 2 = E J pU]r l x l U 1 ^ 1 d z x 

i 



In the ellipsoidal approximation, they are written as 

j, = _±&M WKnp L/n 

1 n + 2 c 

5-n W 2 ,~ 
JL = ; r (A n + n + 2), 

2 3(n + 2) M h 

Using the same method as that adopted in calculating K\, 

R 2 



K 2 = 2j Sxp + 4 + (U 1 - 1 + E XiUy 1 
= ^-W + 2(I un + I U22 ) + 2 J2(Miai + AW 



is rewritten as 



k 3 = -g/^E 



d_ 

dxi 



7U 



i->i 



^2,2 - 



5 



+4xj^— ^TsxDx + 7x 2 -D 2 + x\D 2]2 + ^\D\,\ - Xix 2 D li2 - x x x 2 D 2 ± 
Using the same method as that in calculating J 5 with Eq.( |4.19| ), 



R 2 20 W 

K 3 = T (7W-M 21 )---I, 



3(5 -n) 3(n-4) rj n 

c + 5 - n) 

5 n + 1 K n 



Gathering these results, the PN scalar virial relation is written as 
3 -hK P ^ + l B K 2 pf' = C K ± + D M " 



9 * rc ~ '"92' 

n cr ai <raf 



where 



c = w 



n + 2 



5R 3 



U, 



n 
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5(n - 4) 

n + 1 
135 - 15n 

n+1 



9n + 15 I0rj n 



(5-n) 



60.R 4 



^86 + 35^ -68a 2 ). 
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Using Eqs. (|3~3|) and (|5.15| ), ai = aiN + Sai/c 2 is determined as 



a in 



f>i \"/(n-3)/ £ 



l-n<5M 



v47TQ;2«3|^ / | 
n 



l/(3-n) 



M, 



(n-l)/(n-3) 
N ) 



+ 



Mn 



7T.X 2 

h 



(5.17) 



.3 — n M N ' C{n — 3) a 1N 

where we split M as Mn + 5M/c 2 . The first term of the equation of 5a\ originates from the 
Newtonian property, and the second and third terms (which include B and D) are purely 
the PN origin. Since B < and D>0for0<n<3, the PN effect tends to decrease 8a\ 
for realistic neutron stars, and as a result makes the coordinate radius of the star smaller 
than that in the Newtonian case. 



VI. ANALYSIS OF EQUILIBRIUM SEQUENCES 



An equilibrium sequence of the compressible Darwin ellipsoids is determined fixing M* 
and K. The procedure is as follows. 

(1) using Eqs.( |3.15| ), we calculate the equilibrium sequence in the Newtonian order. Up to 
this stage, a 2 (R) and a 3 (R) are determined, and the results are checked by comparing those 
by LRS @. 

(2) M is determined by the condition M* =constant. From Eq. ( |3.12| ). 



M M. 



+ — + — S-(15gu +4g s 



(6.1) 



Then, we set Mn = M* =constant and 



0-1N I 
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AR 20R 3 
20 



(15g n +4g s ) . 



(6.2) 



(3) after substituting Eqs. (|5.17l ) and ( |6.1| ) into the equations for Q 2 , E\ and J\, we rewrite 
them as 



n 2 



Ei 



a lN 

Ml 

0-1N 



cr I V M* aiN / aiN 



Oin 



2 

PN 



J x = M 3 / 2 a^ N 



<^n H — n 



1 f/35M . 5ai , M, 



+ 



2<2lN 
+ 



4 



M * v 



PN 



aiN 21^ 



Jn H Jpn 

aiN 



(6.3) 
(6.4) 
(6.5) 



. 2M* 2ai N 

and define the total energy and the total angular momentum of the system as E = 1E\ and 
J = 2J\. Then, using ot2{R) and a^^R) obtained at (1), we calculate the PN sequences of 
Q, E, and J as a function of R. 

We repeat these procedures changing the compactness of each star of binary defined as 
C s = M^/^a*, where a* is the radius of the Newtonian spherical star 



M in-i)/{n-3)^ (!^+± K \ n/{3 n) ( 4 < 2 |^|)- (1 - n)/(3 - n) . 



47T 



(6.6) 



For the following, we also define 



3M« 
Anal 



(6.7) 



In figs.l and 2, we show the energy (E /(M 2 /a^)), the angular momentum ( J/ y M% a*), 
and the angular velocity (Q/^/Wp^) as a function of the orbital separation, R^/a* for 
n = 0.5(figs.l) and l(figs.2), respectively. In figs.l, the solid and dotted lines denote the 
results calculated by the ellipsoidal approximation for C s = and 0.02, and the open and 
filled circles denote the results by the numerical calculation ||J] for C s = and 0.02, respec- 
tively. In figs. 2, the solid and dotted lines denote the results calculated by the ellipsoidal 
approximation for C s = and 1/60, and the open and filled circles denote the results by 
the numerical calculation for C s = and 1/60, respectively. These figures indicate that the 
ellipsoidal approximation works fairly well for calculating the angular momentum and the 
angular velocity, but does not very well for calculating the energy: The energy calculated 
by the ellipsoidal approximation systematically deviates from the numerical results. This is 
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because we do not rigidly take into account the change of the density profile of the star due 
to the PN gravity even for R — > oo, although the energy is sensitive to the density configu- 
ration^. e., W or Aq). The angular momentum and the angular velocity are not sensitive to 
the density profile, but to the global quantities such as the mass and the spin as it can be 
found from their equations. However, the deviation of the energy is systematic, so that this 
weak point could be improved if we artificially reduce Aq by a small correction factor which 
will depend on C s and rzfl. 

The ellipsoidal approximation is not quantitatively a sophisticated model for calculating 
the various quantities, in particular the energy, for close binaries of R ~ 3. Nevertheless, 
it is a qualitatively good model: Figs.l and 2 indicate that the coordinate separation at 
the ISCCO becomes small due to the PN effectfl, and this feature agrees with the previous 
numerical results [|J. We show the relation between the angular momentum and the angular 
velocity for n = 0.5 and 1.0 in figs. 3. These indicate that Q at the ISCCO increases with 
the PN effect, and this feature also agrees with the numerical results 

In figs. 4— 8, we show the angular momentum as a function of i?*(a) and the angu- 
lar velocity as a function of the angular momentum(b) for n = 0.25(figs.4), 0.5(figs.5), 

0. 75(figs.6), 1.0(figs.7), 1.5(figs.8) and for C s = 0(solid line), 0.01(dotted line), 0.02(dashed 
line), 0.03(long dashed line), respectively. We can see the feature found in figs.l— 3 again; 

1. e., (1) the coordinate radius at the ISCCO decreases with the increase of C s , and (2) the 
orbital angular velocity at the ISCCO increases with C s . Thus, we conclude that such fea- 

1 For the same C s , we will need a larger correction factor for larger n, and for n = 0, we do not 
need the correction. 

* Since we use the ellipsoidal approximation for the density profile which is not the strict one, and 
also neglect sufficiently high order terms of R~ l in the energy and the angular momentum, the 
location of the energy minimum does not coincide with that of the angular momentum minimum. 
However, the difference is not large like the results in a previous paper |l|]. 
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tures hold irrespective of the EOS if the ISCCO exists at the Newtonian order. These figures 
also indicate that the star of a softer EOS is more susceptible to the PN effect, and as a 
result the orbital radius and the angular velocity at the ISCCO change much. 

Finally, we argue the stability of each star against the collapse due to the influence of 
the GR gravity from the companion star considering the central density p c . From Eqs.( |3.3j ), 
Q5.171 ), and (|6.1| ), we obtain 



p c = Pn 
= Pn 



1 /5M 35ai\ 

1 + " (l 9 + An - 3n» + & = ^ ^ )W 

c^Oin 1 3(3 — n) \ n + 1 



+4i+4i+o( J R- 5 )} 



(6.8) 



where 



47r|^ / |a2a3df N 



^ = , » ( 6 - 9 ) 



and A 3 and A 4 are functions of R. Note that in W, the effect of 0(R~ 3 ) is implicitly included 
because a 2 and a 3 behave as 1 + 0(R~ 3 ). The lowest order terms of _R(the terms of 0(R°)) 
are concerned in the well known GR instability to the radial collapse for each star itself 
|T0| : Since the coefficient of W is always positive for < n < 3, in the case of n — > 3 or of 
sufficiently large W, p c al/c 2 ~ M*/(c 2 a*) becomes larger than unity, which means that the 
star is swallowed inside its Schwarzschild radius to be a black hole due to the GR gravity. 

In Eq. (|6.8|), the terms of 0(i? _1 ) do not appear: As Wiseman points out, this is the 
consequence of the strong equivalence principle |TT] . Thus, the effect due to the GR gravity 
of the companion star appears from the terms of 0(R~ 3 ) which are concerned in the tidal 
effect in the PN order. As Lai pointed out ||12|| , the Newtonian tidal effect tends to stabilize 
each star; i.e., it reduces the central density. Actually, dp^/dR is always positive in the 
present case. Here, we demonstrate that the PN tidal effect also reduces the central density. 
In figs. 9, we show p c /pN as a function of the orbital separation for n = 0.5(a), 1(b), and 
1.5(c). In each figure, the solid, dotted, dashed, and long dashed lines denote the cases for 



23 



C s = 0, 0.01, 0.02 and 0.03, respectively. These figures indicate that the GR effect due 
to the companion star reduces the central density. This seems to mean that it will act to 
stabilize each star of binary. Recently, Wilson et al. calculated the equilibrium state of 



binary neutron stars using the semi-relativistic approximation ||13|| , and they showed that 
the central density of each star increases with the decrease of the orbital separation for a 
critical orbital separation. However, the present analysis suggests that for the corotating 
binary neutron stars, such an effect cannot be explained by the first PN correction of general 
relativity. 

VII. SUMMARY 

Using the ellipsoidal model for the density configuration (the ellipsoidal approximation), 
we have calculated the equilibrium states of corotating binary stars of the stiff polytropic 
EOS in the first PN approximation of general relativity. Comparing the equilibrium se- 
quences with those obtained by numerical calculation [|j], we confirm that the ellipsoidal 
approximation is a qualitatively good model to see the several PN effects in binary neu- 
troOAn stars, and then get the conclusions as follows. 

(1) Irrespective of the polytropic index n, the PN effects make the orbital separation at the 
ISCCO smaller and the orbital angular velocity at the ISCCO larger. 



(2) As in the Newtonian case fll2 |, the PN tidal effect from the companion star reduces the 
central density of each star. Thus, it seems to act the stabilization of each star. 

These results will be helpful to check numerical results in PN and/or full GR numerical 
computations of obtaining the equilibrium states of binary neutron stars. 
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Table I 



Quantities for the polytropic stars of different polytropic indices. 



n 


6 


\6'\ 


h 




On 




Vn 







2/V6 





1 


2/35 


8/35 


6/35 


0.25 


2.59209 


.633717 


.108069 


.903969 


.0561363 


.252922 


.156543 


0.5 


2.75270 


.499997 


.229390 


.814828 


.0538683 


.273988 


.143283 


0.75 


2.93452 


.398094 


.360772 


.731532 


.0507645 


.292471 


.131209 


f.O 


3.14159 


.318310 


.500000 


.653455 


.0471121 


.308844 


.120046 


1.5 


3.65375 


.203301 


.795872 


.511500 


.0389291 


.336503 


.0997846 
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Figure captions 

Figs.l The energy(E/(M^/a^), (a)), the angular momentum( J/ ^ M^a*, (b)), and the an- 
gular velocity (Q / ^/vrp* , (c)) as a function of the orbital separation (i?*/ a*) for n = 0.5. 
The solid and dotted lines denote the results calculated by the ellipsoidal approxima- 
tion for C s = and 0.02, and the open and filled circles denote the results by the 
numerical calculation for C s = and 0.02, respectively. Note that the innermost point 
of the sequence is determined from the condition 2a\ = R. 

Figs. 2 The same as in figs.l, but for n — 1 and C s = 0, 1/60. 

Figs. 3 The relation between the angular momentum and the angular velocity for n = 0.5(a) 
and 1.0(b). The meaning of the lines and marks is the same as that in figs.l or 2. 

Figs. 4 The angular momentum as a function of R*/a* (a) and the angular velocity as 
a function of the angular momentum (b) for n = 0.25, and for C s = 0(solid line), 
0.01(dotted line), 0.02(dashed line), 0.03(long dashed line). 

Figs. 5 The same as figs. 5, but for n = 0.5. 

Figs. 6 The same as figs. 5, but for n = 0.75. 

Figs. 7 The same as figs. 5, but for n — 1.0. 

Figs. 8 The same as figs. 5, but for n — 1.5. 

Figs. 9 Pc/pn as a function of the orbital separation for n = 0.5(a), 1.0(b), and 1.5(c). In 
each figure, the solid, dotted, dashed, and long dashed lines denote the results for 
C s = 0, 0.01, 0.02 and 0.03, respectively. 
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